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ABSTRACT 



g_) I We present new, non-abelian, solutions to the equations of motion 

which describe the collective excitations of a quark-gluon plasma at 



^ . high temperature. These solutions correspond to spatially uniform 

H ; 

^ , color oscillations. 
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The long wavelength excitations of a quark-gluon plasma are collective excitations 
which are described by nonlinear equations generalizing the classical Yang- Mills equa- 
tions in the vacuum. Most studies of such equations have been limited so far to their 
weak field limit, where the modes reduce essentially to abelian-like plasma waves. 
The purpose of this letter is to present new, truly non-abelian, solutions that we have 
obtained recently. At leading order in the gauge coupling g, (we assume 5^ <^ 1 in the 
high temperature deconfined plasma), the collective dynamics is entirely described 
by a set of effective equations for the soft gauge mean fields A^{x) which describe 
the long wavelength (A ~ ^/gT) and low frequency {uj ~ gT) excitations ( T denotes 
the temperature) |l|, 0. (Throughout this work, the greek indices refer to Minkovski 
space, while the latin subscripts are color indices for the adjoint representation of the 
gauge group SU{N)). The equations satisfied by A'^^x) are 

[D'',F^,ix)l=j;i'ixl (1) 

where D^' = d^" + ic/A^(x), (A^ = A^T''), and Ff"^ = [D'',D'']/{ig) = F^^T". The 
induced current j""^ describes the response of the plasma to the color gauge fields A^. 
It is proportional to the fluctuations in the phase-space color densities of quarks and 
gluons. Its expression is|l|, |^ 



jTai^) 



3ulJ — v^Wa{x;v), (2) 



for retarded boundary conditions (A^(x) — i> as Xq — ^ —00). The notations here 
are as follows: uj"^ = (2A^ + Nf)g'^T'^/18 is the plasma frequency, v^ = (1, v), where 
V = k/A; is the velocity of the hard particle with momentum k (A; = |k|), and 
the integral / dVt runs over all the directions of the unit vector v. Furthermore, the 
functions Wa{x\ v) are generally nonlocal and nonlinear functionals of the gauge fields, 
defined by 

/■oo 

Wa{x]v) = I duUab(x,x — vu)v -Ebix — vu) (3) 

Jo 



where E^ = F*° is the chromoelectric field, and U{x, y) is the parallel transporter along 
the straight line 7 joining x and y. The induced current (|^) is covariantly conserved, 
D'^, JT'^{x) = 0. The energy density of an arbitrary gauge field configuration in the 
plasma has been recently computed as[0 

T'^x) = l(E.(a:)-E,(x) + B„(x) • B„(x)) + ^ujI J ^W'Ax;v)W'Ax;v), (4) 

with Bl{x) = —{\/2)t^^^Fl^[x). Remark that T°*'(x) is manifestly positive, and so 
is therefore the excitation energy £{t) = J d^xT^^{t,x.), at any time t. 

The solutions to the non-abelian field equations ([I|) have direct physical relevance: 
they correspond to collective excitations of the high temperature quark-gluon plasma. 
In this letter, we study particular solutions which are uniform in space. (Another 
interesting limit, that of a static field configuration, has been considered recently^, 
with the conclusion that no finite energy solution exists). For convenience, we choose 
the temporal gauge, A°(a;) = 0. Hence, A^{x) = (0, Aa(t)). 

For uniform fields, the functions Wait; v) are simple local functional of the gauge 
potentials 0, 

Wa{t;v) = -v-A,(t), (5) 

and the same holds for the induced current (H), j™'^(t) = — cj^ Aa(t), p^^'^it) = 0, 
and for the energy density (^, 

T'\t) = ^ (^ • ^ + ^>. ■ A.) + ^ r"^ r'^ (A, . A,) (A, . A.) . (6) 



2\ dt dt P V 4 

We have used here the expressions of the field strengths in terms of the vector poten- 
tials, 

Ea(t) = -^, B„(t) = ^r'^Mt) X A,(t). (7) 

The field equations (|^) reduce then to the following equations for the vector 
potentials A^{t) = A^(t)T", 
d'^A' 



dt^ 



+ i^pA' + g' 



A\A^ 



A^ 



0, 



together with a constraint (Gauss's law) 



A' 



dt 



0, (9) 



corresponding to the component yU = of eq. (|I]). These equations are similar to 
the classical Yang-Mills equations in the vacuum, which have been extensively in- 
vestigated already for the case of S't/(2)|^, H, 0, §|. They differ, however, by the 
presence of the thermal mass term ujpA\ which, as we shall see, has a strong effect 
on the dynamics. Note incidentally that, for gauge fields satisfying Gauss's law @, 
the Poynting vector S^ = T*° vanishes 0. 

The constraint (^) is satisfied, in particular, by field configurations of the form 
y4^(t) = Al^h^{t) (no summation over i), with constant Al^ and arbitrary functions 
h^{t). Indeed, for such fields, the three color vectors {A^(t)} and {dAl^/dt} {i = 
1, 2, 3) are parallel in color space. In the rest of this letter we restrict ourselves to 
SU{2), and assume Al^ = 5^. This Ansatz is equivalent to that proposed by Baseyan 
et al.0, up to a trivial global gauge rotation [§]. The functions hi{t) then satisfy 



d'hi 



+ uj'h + g'hihi + hi) = 0, (10) 



dt^ f 

plus two similar equations for /i2 and h^. The associated energy density, 

T°° = ^e((^)' + ^J^.^) +y(^?^2 + /^?/^^+/^>^), (11) 

is an integral of motion and acts as an effective Hamiltonian for the functions hi{t). 
At this point, it is convenient to make a scale transformation and define the dimen- 
sionless variable x = ujpt, as well as the dimensionless functions fi{x) = {g/ujp) hi{t). 
We also assume /s = 0, with no significant loss of generality. We obtain then for /i 
and /2 the coupled nonlinear equations 

/i(x) + [i+{ux)y]Mx) = 0, 

hix) + [l+(/i(x))']/2(x) = 0, (12) 



where the overdots indicate derivatives with respect to x. The energy density (|TT]) is 
then 

4 -, 4 

^r> J- / ^o Ao oO oO oO oO \ ^Ti 



9 ^ ^ ^9 

The Hamiltonian 7i is that of a system of two nonhnearly coupled harmonic oscil- 
lators with coordinates /i and A- Note that the different parameters characteriz- 
ing the initial system (i.e., Up, g, T°°) combine into a single, dimensionless, one, 
6"^ = {g"^ / ujp)T^^ , which measures the total energy of the mechanical system: 7i = 6"^. 
In deriving eq.([l|), we have assumed the gauge fields to be weak {A ^ T) and 
slowly varying {dA ~ (7TA)[|l|. Remembering that ujp ~ gT, one sees that these 
conditions imply |/j(x)| ^ 1 and |/i(a;)| ~ \fi{x)\, so that 9 is, at most, of order unity. 



6^1. These limitations are consistent with the dynamics described by eqs. (|T2D. 
The quadratic terms /^ + /| in the Hamiltonian (|13|), which originates from the 
thermal mass, play an essential role in this respect. Because of them, and of energy 
conservation, a trajectory {fi{x)} cannot leave the bounded doTaain delineated by the 
equipotential lines /i + /| + ff /| = 29^ (Fig.l). Let us assume 9^1. Then, since 
|/j(a;)| < ^/29 (see Fig.l), it follows that |/i(a;)| ^ 1 for any x. Furthermore, the 
quantity (1 + /f (x)), which plays the role of an effective frequency squared for the 
motion in the direction j y^ i, remains of order 1 for any x, so that /j and /j remain 
of the same order of magnitude. Consequently, if the conditions on the functions 
fi mentioned above are satisfied for some xq, then they are valid for any x. The 
situation here is different from the classical, vacuum, case, were the equipotential 
lines are given by hf h\ = {2/g'^)T^^] then, for given energy T^^, the motion could 
extend arbitrary far along the hi or the /i2 axis. Using the language of dynamical 
systems, one may observe that, for the system (0), the origin of the four- dimensional 
phase-space (/i, /i, A, A) is an elliptic fixed point, corresponding to a neutrally stable 
equilibrium [|10|. At T = 0, this same point is marginally unstable. Thus, apart from 
ensuring bounded trajectories {fi{x)}, the quadratic terms f^ in eq. (|T3p also improve 



the stability of the system. 

Let us consider now simple, periodic, solutions to the system (|I2|). 

(a) The simplest motion is one-dimensional along the axes 1 or 2. Assume, e.g., 
/2 = 0. We get then a simple harmonic oscillator equation for /i, 

A(x) + /i(x) = 0. (14) 

The general solution is /i = ai cosa;+a2 sinx, with af+a^ = 26"^. The corresponding 
field configuration, A^(t) = {Cicosujpt + C2smujpt)T^, (with Ci = {ujp/g)ai), and 
A? = A? = 0, describes periodic oscillations (with period Tq = 27i/ujp) in space-color 
direction 1. 

(b) Another one-dimensional, but less trivial, example corresponds to periodic 
solutions with f2{x) = ±/i(x) = ±/(x). These solutions describe in or out of phase 
oscillations of the colors 1 and 2. The function f{x) satisfies the nonlinear equation 

fix) + fix) + fix) = 0, (15) 

for which the integral of motion is 

9' = r + f + l f. (16) 

A particular solution to ([T5| ) is 

fix) = fecn{i29^ + l)'/\x-xo); k), (17) 

where cn(x; k) is the Jacobi elliptic cosine of argument x and modulus fc, and Xq is 
the arbitrary origin of the time. The parameters k and fe are related to 9 by 



1 . 1 vn 



* = 71 1^ - TIFTT I • '^«' 



and 

1/2 



fe = (VWTl -l). (19) 



It can be easily seen that /i = ±/2 = i/g are the coordinates of the intersection 
points between the trajectories and the equipotential hnes in Fig. 1. The solution 
( |T7| ) corresponds to the initial conditions /(xq) = fe and /(xq) = 0. For the in-phase 
oscillations, the associated gauge potentials are A^it) = h(t)T^, v4^(t) = h(t)T^, 
and A^(t) = 0, where h(t) = (ujp/g) fiojpt) is a periodic function, with period 

and K{k) is the complete elliptic integral of modulus k. Since 6 ^ 1, T remains of 
order of Tq = 2TT/ujp. Thus, it is the plasma frequency Up ~ gT which controls the 
time variation of the nonlinear color oscillations (p!?!), for any value of the parameter 
6^1. This is true, in particular, for small oscillations, that is, when the energy 
density T°° ^ 0, so that 6* ~ /^ ^ 1 and /c — >• 0. Then, the periodic solution 
( p!7| ) reduces to a simple harmonic oscillation, and T —>■ %. In contrast, in the zero 
temperature case, as T^^ -^ 0, the corresponding period is diverging[^]. (The T = 
expression for T is obtained from eq. (^) by replacing (29'^ + lY^^cUp -^ (2(7^T°°)^/^ 
and k -^ 1/V2.) 

For this solution, the non-trivial field strengths are i?* = —5^ h, with i = 1, 2, and 
Bl = g h"^ 6^. Accordingly, the three vectors Ei, E2, and B3 are mutually orthogonal. 

An important question concerns the stability of the periodic orbits of the sys- 



tem ([TJ). For the corresponding equations at zero temperature, it has been shown, 
through a numerical analysis, that unstable periodic trajectories exist 0. In particu- 
lar, the in-phase oscillations of two colors (the vacuum analog of i^l7\)) are unstable]^, 
^. In the case of the system (0), the presence of the linear terms /i and /2, which 
originate from the thermal mass in ([TT|), guarantees the stability of the solutions, 
for small oscillations and for most initial conditions. This is verified in the numerical 
analysis of the system (0) performed in Ref. in the physically different context of a 
classical Yang-Mills-Higgs system in the vacuum. There, the mass terms for the gauge 
fields are generated through spontaneous symmetry breaking and are proportional to 



the vacuum expectation value of the scalar field, rj = {(p). The correspondence be- 
tween the two sets of equations is simply given by u^ — » g'^rf /2. The numerical 
experiment in Ref. |0 clearly shows that for ^ ^ 1 (e.g., for 9 ^ 0.2) most of the tra- 
jectories are quasi-periodic. As 9 increases, the stochastic motion strongly develops 
and, beyond a critical value 91 ~ 6.6, it fills the entire permissible range of motion in 
the phase space. As the above value of 9^ stays in the limits of validity of the present 
approach, it would be worth to further investigate the physical signification of this 
transition. Analogous findings, concerning the non-integrability and the stochastic 
behaviour of the dynamical system (p!2[) , are presented in Refs. [|l^ and [Q. 

In conclusion, we have studied here global color oscillations of the high temper- 
ature quark-gluon plasma, by looking at the spatially uniform limit of the effective 
equations of motion derived in Refs. ||l[ . We have shown that the plasma frequency 
ujp = gT{N + Nf/2Y^'^/3 characterizes the inertial properties of the plasma both for 
the abelian-like and for the genuine non-abelian collective motion. The phase-space 
motion is bounded and quasi-periodic for small oscillations and for most initial con- 
ditions. We have obtained explicit nonlinear solutions for SU{2) which describe in or 
out of phase oscillations of two colors (the generalization to the in-phase oscillations 
of the three colors, i.e., hi = h2 = h^ in eq. (|10|), is straightforward]^). Such solutions 



can be imbedded in any larger SU{N) theory in the standard way|T3 
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